Magnetic state of pyrochlore Cd 2 0s 2 07 emerging from strong competition of ligand 
distortions and longer-range crystalline anisotropy 
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By many-body quantum-chemical calculations, we investigate the role of two structural effects - lo- 
cal ligand distortions and the anisotropic Cd-ion coordination - on the magnetic state of Cd20s207, 
a spin S = 3/2 pyrochlore. We find that these effects strongly compete, rendering the magnetic in- 
teractions and ordering crucially depend on these geometrical features. Without trigonal distortions 
a large easy-plane magnetic anisotropy develops. Their presence, however, reverses the sign of the 
zero-field splitting and causes a large easy-axis anisotropy (Z?~ —6.8 meV), which in conjunction 
with the antiferromagnetic exchange interaction ( J ~ 6.4 meV) stabilizes an all-in/all-out magnetic 
order. The competition uncovered here is a generic feature of pyrochlore magnets. 



Introduction. Oxide compounds of the 5d elements 
are at the heart of intensive experimental and theoreti- 
cal investigations in condensed-matter physics and ma- 
terials science. The few different structural varieties dis- 
playing Ir ions in octahedral coordination and tetravalent 
5d 5 valence states, in particular, have set the playground 
for new insights into the interplay between strong spin- 
orbit interactions and electron correlation effects [1-3]. 
Equally intriguing but less investigated are pentavalent 
Os 5d 3 oxides such as the pyrochlore Cd20s207 [4, 5] 
and the perovskite NaOsC>3 [6, 7]. Just as for the iri- 
dates, major open questions with respect to the osmatcs 
are the mechanism of the metal to insulator transition 
(MIT) and the nature of the low-temperature antifer- 
romagnetic (AF) configuration. While initially a Slater 
type picture has been proposed for the MITs in both 
Cd 2 0s 2 07 [4] and NaOs03 [6, 7], alternative scenarios 
such as a Lifshitz transition in the presence of sizeable 
electron-electron interactions have been suggested as well 
very recently [8, 9]. 

So far theoretical investigations of the electronic struc- 
ture of the Os 5d 3 oxides have been carried out only 
as mean-field ground-state (GS) calculations within the 
local-density approximation (LDA) [10] or the LDA+C7 
model [6, 8]. Here, we report results of many-body 
quantum-chemical calculations for Cd20s207. We de- 
scribe the local Os d 3 multiplet structure, the precise 
mechanism of 2nd-order spin-orbit coupling (SOC) and 
zero-field splitting (ZFS), and determine the parame- 
ters of the effective spin Hamiltonian, i.e., the single- ion 
anisotropy (SIA), nearest-neighbor (NN) Hcisenberg ex- 
change as well as the Dzyaloshinskii-Moriya (DM) in- 
teractions. The results indicate that basic electronic- 
structure parameters such as the SIA crucially depend 
on geometrical details concerning both the O and adja- 
cent Cd ions. In particular, the anisotropic Cd-ion coor- 
dination lifts the degeneracy of the Os t2 g levels even in 
the absence of O-ion trigonal distortions and yields large 
ZFSs and easy-plane magnetic anisotropy. A configu- 



ration with trigonally distorted O octahedra is however 
energetically more favorable, changes the sign of the ZFS, 
and thus gives rise to easy-axis anisotropy at the Os sites, 
which in conjunction with the AF NN exchange stabilizes 
the so-called all-in/all-out spin order. The competition 
between the two structural effects is a generic feature in 
227 pyrochlores and opens new perspectives on the basic 
magnetism in these materials. 

Os d 3 electron configuration, single-ion physics. To 
investigate in detail the Os d-level electronic struc- 
ture, multiconfiguration self-consistent-field (MCSCF) 
and multircfcrcncc configuration-interaction (MRCI) cal- 
culations [11] were first performed on embedded clusters 
made of one reference Os06 octahedron, six adjacent Cd 
sites, and six NN Os06 octahedra. The farther solid- 
state environment was modeled as an one-electron ef- 
fective potential which in a ionic picture reproduces the 
Madclung field in the cluster region. All calculations were 
performed with the molpro quantum-chemical software 
[12]. We employed energy-consistent relativistic pseu- 
dopotcntials for Os [13] and Cd [14] and Gaussian-type 
valence basis functions from the MOLPRO library (see 
Supplementary Material [15]). 

From simple considerations based on crystal-field (CF) 




FIG. 1. Sketch of (i) the compressive trigonal distortion in 
Cd20s207 (red arrows) involving the O ligands (in beige) 
around each Os site (blue) and (ii) the hexagonal Cd-ion 
(green) coordination. The clusters used in the calculation 
of the ZFSs also include the six NN octahedra, see text. 
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TABLE I. MRCI and MRCI+SOC relative energies (eV) for 
the Os 5+ 5d 3 multiplet structure in CC12OS2O7. Since cubic 
symmetry is lifted, the T states are split even without SOC. 
Each MRCI+SOC value stands for a spin-orbit doublet; for 
the 4 T states only the lowest and highest components are 
shown. 
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theory, the orbital GS for a d 3 ion in octahedral coor- 
dination is singlet. Thus, neglecting SOC, the lowest 
electron configuration should be 4 A2(t\ g ). This is indeed 
found in the quantum-chemical calculations if SOC is 
not included. The components of the spin-quartet Ai 
GS can interact, however, via SOC with higher lying 
T2 terms and for noncubic axial systems split into two 
Kramers doublets, mj = ±1/2 and wig = ±3/2, respec- 
tively [16, 17]. The T2 states are themselves split in non- 
cubic CFs. In the recent work of Matsuda et al. [5], for 
instance, the interaction of the 4 A2 and lowest 4 T2 states 
was analyzed. The ab initio quantum-chemical calcula- 
tions yield, nevertheless, a rather large ^A^-^T^, excita- 
tion energy, 5.4 eV by MCSCF calculations including in 
the active orbital space all 5d functions at the central Os 
site and the t2 9 orbitals of the six Os NNs, which is not 
surprising in fact for extended 5d orbitals that feel very 
effectively the O 2p charge distribution. Additionally, the 
ab initio investigation actually shows that the excited 2 E 
and 2 T t^g doublets occur at much lower energies, 1 to 2.5 
eV (see Table I) , and the ZFS of the 4 A2 levels is mainly 
related to the interaction with those states. To make the 
MRCI calculations feasible [18] and the analysis of the 
results limpid, we further replaced the six Os 5+ d 3 NNs 
by closed-shell Ta 5+ d° species. We checked, for example, 
that the bd charge distribution and the t2 g -e g splitting at 
the central Os site are not affected by this substitution. 
In the MRCI computations, we included on top of the 
MCSCF wave functions all single and double excitations 
from the Os 5d and O 2p orbitals at the central octa- 
hedron. The reference multiconfiguration active space is 
given by five 5d orbitals at the central Os site and three 
electrons and the results are listed in Table I. 

To extract the ZFS, we compute the spin-orbit interac- 
tion matrix for the lowest four quartet and three doublet 
spin states, see Table I. The MRCI+SOC data of Ta- 
ble I show that the ZFS of the 4 A2 GS is rather large, 
13.5 meV. This gives a first estimate of the splitting pa- 
rameter \D\ ~6.75 meV, one to two orders of magnitude 



larger than for 3d ions in standard coordination [16, 17]. 
However, for explicitly deriving the full SI A tensor D, 
we further employ the effective Hamiltonian methodology 
described in Rcf. [19]. Here, SOCs and the mixing of the 
4 A2 components with higher lying CF excited states are 
treated as small perturbations and the spin-orbit wave 
functions related to the high-spin t^g configuration are 
projected onto the space spanned by the 4 ^2 \S,Ms) 
states. Using the orthonormalized projections tyf. of the 
low-lying A2 quartet wave functions and the correspond- 
ing eigenvalues we then construct the effective Hamil- 
tonian n cS = Y.k E k\^k){^k\ [19]. The orthonormal- 

ization is done using the formalism introduced by des 

_ A 

Cloizcaux [20], i.e., \^k) = U ov 2 l^fc), where U ov is the 
overlap matrix and \^k) are the projections of the wave 
functions onto the model space. A one-to-one correspon- 
dence can be now drawn between the matrix elements of 
'HeB and the model Hamiltonian for the anisotropic sin- 
gle site problem, 7^ mo d = S • D ■ S, to extract the ZFS 
tensor (for details, see Supplementary Material [15]). Di- 
agonalizing the D tensor we finally obtain on the basis of 
the ab initio quantum-chemical data an axial parameter 
D = —6.77 meV, a rhombic parameter E = 0.0, and an 
univocally defined easy axis of magnetization along the 
[111] direction, perpendicular to the plane of Cd NNs. 

To gain deeper insight into the origin of the strong ax- 
ial anisotropy, we performed additional calculations for 
an ideal pyrochlorc structure with no trigonal distortions 
of the O octahedra. It turns out that the 227 pyrochlore 
lattice is fully defined by just three parameters: the space 
group, the cubic lattice constant a, and the coordinate x 
of the O at the 48/ site, denoted 0(1) in Rcf. [4]. For 
o;=:eo = 5/16, the oxygen cage around each Os site forms 
an undistorted, regular octahedron [21]. In Cd20s 2 07, 
however, x is slightly larger than xq [4], which translates 
into a compressive trigonal distortion of the OsOe octa- 
hedra. Interestingly, for undistorted O octahedra (i.e., 
x = xq) we still found E = and the same orientation 
of the magnetic z axis. The axial parameter, however, 
changes its sign to £> = 2.50 meV. A positive D indicates 
that the z axis is now a hard axis. The basic mecha- 
nism responsible for this change of sign is apparent when 
comparing the Os 5c? 3 multiplet structure for the actual 
lattice (Table I) with that for a idealized crystal with no 
trigonal distortions (Table II). In particular, the order of 
the split components of the 2 T\ states, i.e., one singly- 
degenerate and one doubly-degenerate term, changes in 
Table II as compared to Table I, which modifies the in- 
teraction of the A2 terms with the excited states. 

We find in the MRCI calculations a linear relation be- 
tween D and the distortion angle A#, see Fig. 2. The 
effect of a trigonal field alone on the energy levels and 
on the ZFS of d 3 ions has been previously investigated 
in the framework of a simplified CF model and perturba- 
tion theory [17, 22]. A linear relation is also established 
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TABLE II. MRCI and MRCI+SOC relative energies (eV) for 
the Os 5+ 5d 3 multiplet structure in a hypothetical crystal 
with no trigonal distortions. Each MRCI+SOC value stands 
for a spin-orbit doublet; for the 4 T states only the lowest and 
highest components are shown. 
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in alum salts [22]. In contrast to the model assumed in 
Ref. [22], however, we find that the D(A8) line does not 
go through origin. This is due to the noncubic field gen- 
erated by the Cd NNs, which form a highly anisotropic 
planar structure arround each Os site. Similar competing 
stuctural effects have been evidenced in Ni n -Y m molec- 
ular complexes [23]. We also plot in Fig. 2 the total- 
energy landscape as function of the distortion angle A9. 
It is seen that at the MRCI level the minimum of the 
parabola corresponds to a trigonal compressive distor- 
tion A9 = 1.7° for fixed Os-0 bond length, close to the 
experimental value of 1.9° [4]. 

Superexchange interactions. To investigate the NN 
magnetic couplings, we designed ten-octahedra embed- 
ded clusters. The analysis of the intersite interactions is 
here carried out for only two magnetically active Os d 3 
ions (see Fig. 3), while for simplicity the eight Os 5+ NNs 
were modeled as closed-shell Ta 5+ d° ions. Multicon- 
figuration wave functions were first generated by state- 
averaged MCSCF optimizations for the lowest singlet, 
triplet, quintet, and septet states in the two-site prob- 
lem. Those configuration state functions give rise in the 
spin-orbit calculations to sixteen spin-orbit states. For 
each spin multiplicity the MCSCF active space is defined 
by the set of six Os t2 g orbitals accommodating a to- 
tal number of six electrons. We then further accounted 
for single and double excitations from the Os t2 g and 
bridging-0 2p orbitals on top of the MCSCF reference 
wave functions. Such MRCI calculations yield magnetic 
coupling constants in very good agreement with exper- 
imental data in the 5d 5 iridatc Sr2ix04 [24]. Similar 
strategies of explicitly dealing only with selected groups 
of localized ligand orbitals were adopted in earlier studies 
on 3d compounds [25]. 

From the calculations without SOCs, we see that the 
energy splittings between the different spin states fol- 
low the sequence J, 3 J, 6J and can be fitted to an AF 
isotropic exchange model Haf = JSi-&2, with a mean 
squared error of 0.6 meV. Adding a biquadratic term 
K(Si ■ S2) 2 improves the accuracy of the fit to 0.01 meV 



and yields a first MRCI estimate of J = 6.42 meV. Thus 
J is of the same order of magnitude as the SIA while 
A' = 0.07 meV. 

Antisymmetric exchange. Anisotropic, non- 

Heisenberg terms can be further obtained from 
spin-orbit calculations within the manifold defined 
by the lowest singlet, triplet, quintet, and septet spin 
states for two Os d 3 sites. As expected, including SOCs 
does not bring significant corrections to the effective 
superexchange: J is 6.43 meV by MRCI+SOC calcu- 
lations. The spin-orbit interactions, however, lift the 
degeneracy of the initial spin multiplcts and brings in 
antisymmetric contributions. To describe the latter, we 
adopt a two-site model Hamiltonian 7^ lod containing, 
in addition to the Heisenberg and biquadratic terms, 
a DM antisymmetric exchange contribution d-Si x S2, 
where d is a pseudovector with components d x , d y , d z . 
This model is then compared with the ab initio 16 x 16 
effective matrix TL' cS derived from the MRCI+SOC 
calculations [15]. We found perfect one-to-one corre- 
spondence between the two sets of matrix elements 
and could extract the numerical parameters J = 6.43, 
K = 0.07, and d = (1.17,-1.17,0) meV. According to 
these results, the DM vector is oriented along the (110) 
axes (see Fig. 3 and Ref. [26]) and its norm is |d| = 1.65 
meV. The sense of the DM vector for a given Os-Os link, 
however, cannot be determined from our ab intio data. 

Discussion. Whereas the absolute value of the MRCI 
DM parameter is similar to estimates derived for certain 
values of the on-site Coulomb repulsion U by LDA+J7 
calculations [8], the superexchange J and the SIA D 
turn out to be two (J) to seven (D) times smaller in 
our quantum-chemical study. It is known that quantum- 
chemical calculations based on either configuration- 
interaction techniques or 2nd-order perturbation theory 
are able to reproduce the experimentally derived J cou- 
pling constants with an accuracy better than ±20% in 
insulating 3d-metal oxides, see, e.g., Refs. [27] and [25]. 
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FIG. 2. (a) Single-ion anisotropy parameter D and (b) total 
energy of the cluster A_B to t, as functions of the trigonal dis- 
tortion angle A8, see text. MRCI+SOC results for clusters 
as that depicted in Fig. 1. No trigonal distortions are applied 
within the embedding, i.e., x' = x'o = 5/16 for the surrounding 
lattice. 
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For more extended 5d electronic states in 5d-metal ox- 
ides, the Hubbard U is neither small nor large compared 
to the bandwidth, as corroborated from the MITs [4, 7] 
and from constrained LDA calculations [6, 8, 28]. Elec- 
tron itineracy may therefore pose additional technical 
problems and limit the applicability of a finite cluster 
approach. It has been shown, however, that accurate J's 
[24] and d-d excitation energies [29] can be indeed com- 
puted for 5c? oxides such as the iridates. Those findings 
in [24] and [29] indicate that approaching the essential 
physics in the 5d oxides from a more localized perspec- 
tive and within a finite- fragment framework certainly is a 
good starting point. Additional potentially problematic 
points concern the size of the clusters, as discussed, e.g., 
in Ref. [30]. The clusters we employ here are, however, 
large enough to ensure an accurate description of both 
the charge distribution of NN octahedra around the "ac- 
tive" d sites and the tails at those neighboring sites of 
Wannier-like orbitals in the active region. As concerns 
the SIA at d-metal sites in solids, extensive quantum- 
chemical investigations are missing. The present inves- 
tigation therefore constitutes a step toward filling this 
gap. 

The pure AF Heisenberg model on the pyrochlore lat- 
tice has an extensive number of classical GSs which pre- 
vent any Neel type ordering down to zero temperature 
[31, 32]. The easy-axis anisotropy lifts this macroscopic 
degeneracy and selects the so-called all- in/all-out (or it's 
time reversed) state, whereby all four spins sharing a 
given tetrahedron point toward or away from its center. 
This GS is selected irrespective of the relative strength 
of the anisotropy and AF exchange energy [33]. Now, 
since the macroscopic GS degeneracy is already lifted by 
the SIA, including the smaller DM interactions is not ex- 
pected to drastically modify the physics of the system. 
In fact, if the sign of the d vector is that of the direct 




FIG. 3. Network of corner-shared OS4 tetrahedra (blue color) 
in Cd20s207. O ligands and NN Cd ions around two adja- 
cent Os sites are shown in beige and green. The two possible 
orientations of the DM vector [26] for two adjacent Os ions is 
also shown, along with their respective easy axes. 



type (as defined by Elhajal et al. [26]), the DM inter- 
action also favors the all-in/all-out state and thus does 
not compete with the SIA. The competition arises for the 
opposite sign of the d vector, labelled as indirect type in 
[26], where the DM couplings favor a continuous family 
of coplanar and non-coplanar GSs [26]. 

Conclusions. While trigonal distortions of the ligand 
cage are believed to open the door to topologically insu- 
lating states in j = 1/2 227 iridate pyrochlores [2, 28, 34], 
we here show that the lattice degrees of freedom can fun- 
damentally modify the nature of the magnetic interac- 
tions and ground states in 5d 227 compounds with larger 
magnetic moments. For this we employ ab initio multiref- 
erence configuration-interaction methods, the results of 
which are further mapped onto model Hamiltonians in- 
cluding both isotropic and anisotropic terms. To extract 
the DM couplings, we have expanded the formalism pro- 
posed earlier for S —1/2 cupratcs [35]. For the experi- 
mentally determined crystal structure of Cd20s207 [4], 
we find an AF superexchange coupling J = 6.43 meV and 
easy-axis anisotropy of the same magnitude, D = —6.77 
meV, while the DM parameter is 3 to 4 times smaller. 
The dominant magnetic interactions, J and D, give rise 
to the all- in/all-out AF order [26, 33]. Most remarkably, 
the hexagonal Cd-ion coordination induces electrostatic 
fields that compete with the compressive trigonal distor- 
tion of the O octahedra that is present in the system. If 
the latter is removed the SIA turns positive, i.e., D re- 
verses sign and the system develops easy-plane magnetic 
anisotropy. This obviously implies qualitatively different 
magnetic properties [36] . Our findings thus indicate that 
there is in 227 pyrochlores a transition from unfrustrated 
to highly frustrated magnetism controlled by the amount 
of trigonal distortion. The latter can in principle be var- 
ied by changing either the A-site or B-site ions in the 
A2-B2O7 compounds, replacing for instance Cd(Os) by 
Hg(Re) [21], and is in addition tunable by applying pres- 
sure. For a correct understanding of such magnetic in- 
teractions, ordering tendencies, and transitions the subtle 
interplay between oxygen- and A-ion anisotropic electro- 
static effects is essential. 
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Supplementary material 
Single-ion problem 

To determine the Os d-level electronic structure we em- 
ployed seven-octahedra clusters (see Fig. SI) embedded 
in arrays of point charges fitted to reproduce the crys- 
tal Madelung field in the cluster region. Valence basis 
sets of quadruple-zeta quality [SI] from the MOLPRO [S2] 
library were used for the central Os ion [S3] and triple- 
zeta basis functions were applied for the ligands [S4] of 
the central octahedron and the nearest-neighbor (NN) 5d 
sites [S3]. For the central Os ion we also employed two 
polarization / functions [S3]. For farther ligands in our 
clusters we applied minimal atomic-natural-orbital basis 
sets [S5]. All occupied shells of the NN Cd 2+ ions were 
incorporated in the large-core pseudopotentials and each 
of the Cd 5s orbitals was described by a single contracted 
Gaussian function [S6]. We accounted for all single and 
double excitations from the Os 5d and O 2p orbitals 
at the central OsOg octahedron in the multireference 
configuration-interaction (MRCI) treatment. To sepa- 
rate the O 2p valence orbitals into two different groups, 
i.e., at sites of the central octahedron and at NN octahe- 
dra, we used the orbital localization module available in 

MOLPRO. 

We mapped our ab initio data onto a anisotropic single- 
site magnetic Hamiltonian. The interaction matrix in 
the absence of magnetic fields is shown in Table SI and 
the effective Hamiltonian constructed on the basis of 
the MRCI+SOC data is given in Table SIT The onc-to- 




FIG. SI. The cluster used for the ZFS calculations. Os, O, 
and Cd ions are shown in blue, beige, and green, respectively. 



one correspondence between analogue terms in Tables SI 
and SII gives rise to a system of equations that can be 
easily solved to obtain the single-ion anisotropy tensor 
in the original (crystallographic) coordinate frame. The 
later reads, in meV, 

/ D xx D xy D xz \ ( 1.806 -2.258 -2.258 \ 
D= D xy D yy D yz = -2.258 1.806 -2.258 . 
V D xz D yz D zz J \ -2.258 -2.258 1.806 / 

The diagonalization of D yields the transformation ma- 
trix for rotating the initial coordinates to the magnetic 
framework and also the ZFS parameters, i.e., the axial 
D and the rhombic E components: 

(Dxx \ / 4.064 \ 
Dma g = D YY = 4.064 , 
\ £>zz / V -2.710 / 

D = D ZZ - i(£>xx + Dyy) = -6.774 meV, 
E = ^(D xx - Dyy) = 2 • 10~ 5 meV « 0. 

Two-site problem 

Calculations for the inter-site magnetic couplings were 
performed on embedded ten-octahedra clusters with two 
active Os 5+ 5d 3 sites, as shown in Fig. S2. The basis 
functions used here were as described above. The only 
exception was the O ligand bridging the two magnetically 
active Os sites, for which we employed quintuple-zeta 
valence basis sets and four polarization d functions [S4]. 
The eight adjacent 5c? sites were modeled as Ta 5+ d° 
species. Relative energies for the lowest singlet, triplet, 
quintet, and septet are listed in Table SIII. 

The spin-orbit MRCI calculations and the subsequent 
analysis were performed in this case within the manifold 
defined by a S = 3/2 electron configuration at each ac- 
tive Os site. In other words, for this 4 A2® 4 A2 manifold, 
we do not include second-order couplings through higher 
lying crystal-field excited states. As a result, we describe 
only the diagonal interactions and the first-order con- 
tribution to the Dzyaloshinskii-Moriya (DM) anisotropic 
term. For 3d oxides, it has been shown, nevertheless, 
that the diagonal couplings are the most important, sec 
Refs. [S7, S8]. For the same reason, the above treatment 
does not include the single-ion anisotropies, which we 
have obtained from the single-site calculations described 
above. Only 2p orbitals at the bridging O site were in- 
cluded in the MRCI treatment, as in earlier CI calcula- 
tions on 3d transition-metal compounds [S9]. 

We found that the ab initio results can be perfectly 
mapped onto a model Hamiltonian with bilinear and bi- 
quadratic isotropic terms plus antisymmetric exchange, 
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TABLE SI. Model interaction matrix for ri mod = S ■ D • S and S = §. 



Hmod | 2 > 2 ) _ I 2 ' 2 ) I 2 ' 2 ) I 2 ' 2 ) 

!§,-§| z( D ** + D w) + Z D ** -V3(D xz +iD yz ) ^(D xx -D yy + 2iD xy ) 

;f,-|| -V3(D XZ - iD yz ) 1(D XX + D yv ) + \D tt ^-(D xx ~D yy + 2iD xy ) 

(§,±| ^(D xx - D yy - 2iD xy ) + A/sHi-D** + «A/~) 

(|,|| ^(-D M D WB 2//), , : V3(D XZ - iD yz ) \{D XX + D yy ) + \D Z * 




FIG. S2. The cluster used for computing the intersite mag- 
netic couplings. Os, O, and Cd sites are shown in blue, beige, 
and green, respectively. The NN Os 5+ ions around the two 
central octahedra are modeled as Ta 5+ d° species. 

as presented in Table SIV. The effective Hamiltonian 
constructed using the MRCI+SOC data is shown in Ta- 
ble SV. Small non-diagonal elements on the order of 0.01 
meV are neglected in our analysis. One-to-one corre- 
spondence between the elements in Tables SIV and SV 
confirms that our choice of H'mod is appropriate and al- 
lows to extract all the numerical parameters of the model. 



One should notice that only the differences between the 
diagonal elements can be used in the analysis because the 
trace of an effective Hamiltonian is in general arbitrary. 
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TABLE SII. Single-site effective interaction matrix obtained by MRCI+SOC calculations on seven-octahedra cluster (meV). 



HcS 

31 



(I. 



'2 II 
\ 2 ' 2 I 



|3 _3\ 
I 2' 21 

6.7738 



11,-1) 
3.9110 + 3.9110i 



3.9110 - 3.9110i 
3.9106i 
0.0000 



6.7741 
0.0000 
3.9106i 



I- -) 

I 2 ' 2 / 

-3.9106i 
0.0000 
6.7738 
-3.9110 + 3.9110i 



I- -) 
0.0000 

-3.9106i 
-3.9110 - 3.9110i 
6.7738 



TABLE SIII. Relative energies (meV) of the singlet, triplet, 
quintet, and septet states for two adjacent Os d 3 sites. In 
each case, the singlet is taken as reference. 

MCSCF MRCI MRCI+SOC 



5*tot = 
5tot = 1 



Stot — 2 



<Stot — 3 



0.00 
2.51 



7.75 



16.24 



0.00 
5.98 



18.38 



37.93 



0.00 
6.04 
6.04 
6.44 
18.66 
18.66 
19.01 
19.01 
19.12 
38.51 
38.51 
38.81 
38.81 
38.99 
38.99 
39.05 



TABLE SIV. Model interaction matrix in coupled basis \S tot , Mg ot ) for H' mod = JSi ■ S 2 + K(Si ■ S 2 ) 2 + d(Si x S 2 ) and Si, 2 = § 
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TABLE SV. Effective interaction matrix in coupled basis \S tot , M l s ot } , obtained from two-sites MRCI+SOC calculations for Cd 2 0s 2 7 (meV) 





|3, —3) 


|3,-2) 


13,-1} 


|3,0} 


|3,1) 


13,2} 


|3,3) 






(3, -3| 


38.525 


0.000 


0.016i 


0.000 


-0.011 


0.000 


0.000 






(3, -2| 


0.000 


38.505 


0.000 


0.000 


0.000 


0.000 


0.000 






<3,-l| 


-0.016i 


0.000 


38.511 


0.000 


O.OlOi 


0.000 


-0.011 






<3,0| 


0.000 


0.000 


0.000 


38.507 


0.000 


0.000 


0.000 






(3,1| 


-0.011 


0.000 


-0.01? 


0.000 


38.511 


0.000 


0.016; 






<3,2| 


0.000 


0.000 


0.000 


0.000 


0.000 


38.505 


0.000 






<3,3| 


0.000 


0.000 


-0.011 


0.000 


-0.016i 


0.000 


38.525 






(2, -2| 


-2.120 + 2.120i 


0.000 


-0.550 - 0.550i 


0.000 


0.000 


0.000 


0.000 






<2,-l| 


0.000 


-1.727 + 1.727i 


0.000 


-0.947 - 0.947i 


0.000 


0.000 


0.000 






(2,0 


0.000 


0.000 


-1.339 + 1.339i 


0.000 


-1.339 - 1.339; 


0.000 


0.000 






<2,1| 


0.000 


0.000 


0.000 


-0.947 + 0.947i 


0.000 


-1.727 - 1.727; 


0.000 






(2 2| 


000 


000 


000 


000 


— O 550 -1- O 5507 


000 


-9 19(1 - 9 1 90; 






<i,-i| 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 






<1,0| 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 






(1,11 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 






(0,0| 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 






tt'cff 


12,-2) 


12.-1} 


12,0} 


12,1} 


12,2} 


11,-1) 


11,0} 


11,1} 


|0,0> 


(3, -3| 


-2.120 - 2.1202 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


<3, -2 


0.000 


-1.727 - 1.727i 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


<3,-l| 


-0.550 + 0.550i 


0.000 


-1.339 - 1.339i 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


(3, 0| 


0.000 


-0.947 + 0.947i 


0.000 


-0.947 - 0.947i 


0.000 


0.000 


0.000 


0.000 


0.000 


<3,1| 


0.000 


0.000 


-1.339 + 1.339i 


0.000 


-0.550 - 0.550; 


0.000 


0.000 


0.000 


0.000 


<3, 2| 


0.000 


0.000 


0.000 


-1.727 + 1.727; 


0.000 


0.000 


0.000 


0.000 


0.000 


(3,3 


0.000 


0.000 


0.000 


0.000 


-2.120 + 2.120'i 


0.000 


0.000 


0.000 


0.000 


(2, -2| 


18.954 


0.000 


-0.019i 


0.000 


0.000 


1.827 + 1.827; 


0.000 


0.000 


0.000 


<2,-l| 


0.000 


18.977 


0.000 


-0.023; 


0.000 


0.000 


1.293 + 1.293i 


0.000 


0.000 


<2, 0| 


0.019i 


0.000 


18.985 


0.000 


-0.019i 


0.745 - 0.745; 


0.000 


0.745 + 0.745i 


0.000 


<2,1| 


0.000 


0.023; 


0.000 


18.977 


0.000 


0.000 


1.293 - 1.293; 


0.000 


0.000 


(2, 2| 


0.000 


0.000 


0.019; 


0.000 


18.954 


0.000 


0.000 


1.827 - 1.827i 


0.000 


<1.-1| 


1.827 - 1.827i 


0.000 


0.745 + 0.745i 


0.000 


0.000 


6.528 


0.000 


0.000 


0.936 + 0.936i 


<1,0| 


0.000 


1.293 - 1.293; 


0.000 


1.293 + 1.293i 


0.000 


0.000 


6.528 


0.000 


0.000 


(Ml 


0.000 


0.000 


0.745 - 0.745i 


0.000 


1.827 + 1.827; 


0.000 


0.000 


6.528 


0.936 - 0.936i 


<0, 0| 


0.000 


0.000 


0.000 


0.000 


0.000 


0.936 - 0.936; 


0.000 


0.936 + 0.936i 


0.579 



